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^SJ . We report a detailed account of the phase diagram of a recently introduced model for non- 

equilibrium wetting in 1 + 1 dimensions [Phys. Rev. Lett. 79, 2710 (1997)]. A mean field approxi- 
mation is shown to reproduce the main features of the phase diagram, while providing indications 
, for the behaviour of the wetting transition in higher dimensions. The mean field phase diagram is 

^ ' found to exhibit an extra transition line which does not exist in 1 + 1 dimensions. The line sepa- 

' rates a phase in which the interface height distribution decays exponentially at large heights, from 

a superexponentially decaying phase. Implications to wetting in dimensions higher than 1 + 1 are 
discussed. 
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I. INTRODUCTION 



Wetting phenomena occur in a large variety of experiments, where a planar substrate is exposed to a gas phase 
"j^ [ under thermal equilibrium conditions. Generally, 'wetting' refers to a situation where a bulk phase (a) in contact 
with a substrate coexists with a layer of a different phase (/3) which is preferentially attracted to the surface of the 
i' ' substrate. By changing physical parameters such as temperature and chemical potential, the system may undergo a 
, wetting transition from a non-wet phase, where the thickness of the layer stays finite, to a wet phase, where the layer 
^ ■ becomes macroscopic. 

o ■ 

^ ' The phase diagram associated with the surface layer could be rather complex exhibiting a variety of surface phase 
transitions, prewetting phenomena and multicritical behavior 1, 2]. For example, by increasing the temperature T 
^ ' while moving along the a-(3 coexistence curve, a wetting transition may take place at a temperature Tw, beyond 
^ . which the thickness of the layer becomes infinite. Usually this transition is first order, although in certain models 
(~«^ ' the transition is continuous, and is then referred to as continuous wetting. On the other hand, when the chemical 
lO ', potential difference between the two phases is varied, moving towards the coexistence curve at T > Tw, a different 
CO . type of transition takes place in which the thickness of the layer diverges. This phenomenon is known as complete 
^'^^ ' wetting. 

I In many experimental situations it is reasonable to assume that a wetting stationary layer is in thermal equilibrium. 
' In fact, methods of equilibrium statistical mechanics turned out to be very successful in a large variety of theoretical 
and experimental studies (for a review, see Ref. 0). Within this approach, a wetting transition is usually considered 
as the unbinding of an interface from a wall. The interface configuration is described by a function /i(x) which gives 
the height of the interface at point x on the substrate. One then introduces an effective Hamiltonian of the form ^ 



o 



(1) 



where cr is the effective surface tension of the a-[3 interface, V{h) is a potential accounting for the interaction between 
the wall and the interface, and d — 1 is the interface dimension. In the non-wet phase the potential V contains an 
attractive component which binds the interface to the wall. Assuming thermal equilibrium, the probability of finding 
5— i ■ the interface in a certain configuration is then given by the canonical distribution 

P[h] -exp(-/3H[/i]). (2) 

As the parameters describing the system are varied, the attractive component of the potential may become weaker so 
that it is no longer able to bind the interface, leading to a wetting transition. 

In order to study wetting phenomena under thermal equilibrium conditions one usually introduces a stochastic 
Langevin equation corresponding to the effective Hamiltonian This Langevin dynamics should reproduce the 
equilibrium distribution ^ in the limit i — > oo. Since many different dynamical rules may approach the same 
stationary state, this condition does not fully determine the form of the Langevin equation. However, assuming 
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short-range interactions and keeping only the most relevant terms in the renormalization group sense, one is led to 
the Edwards- Wilkinson equation with a potential 0| 

where C(x,i) is a zero-average Gaussian noise field with a variance 

(C(x, t)C(x', t')) = 2r5'^-i(x - ^')d{t ~ t') , (4) 

and a noise amplitude T = ksT. This Langevin equation has the same symmetry properties as the Hamiltonian l(T]l. 
namely, it is invariant under translations, rotations, and reflections in space. Apart from the potential term, the 
equation is also invariant under shifts h ^ h + a and reflections h — > —h. Moreover, it can be shown that this type of 
Langevin dynamics obeys detailed balance and relaxes towards the equilibrium distribution (jSJ in the bound phase. 

Wetting phenomena may also take place in many systems under non-equilibrium conditions. For example, in 
growth processes such as molecular beam epitaxy or others, a layer is grown on a substrate, whose properties depend 
on the growth conditions. By varying these conditions one expects wetting phenomena to take place. Here, unlike 
the equilibrium case, the dynamics does not obey detailed balance, leading to a rather different class of wetting 
phenomena. 

The simplest way to study nonequilibrium wetting on the level of the Langevin equation is to introduce a nonlinear 
term in Eq. leading to a Kardar-Parisi-Zhang (KPZ) equation with a potential Q 

It is important to note that this nonlinear term is a relevant perturbation of the underlying field theory, i.e., even if A 
is very small, it will be amplified under renormalization group transformations, driving the system away from thermal 
equilibrium. 

Recently, a simple solid-on-solid (SOS) model for non-equilibrium wetting in 1 4- 1 dimensions was introduced 0,|3|. 
The model is controlled by an adsorption rate q and a desorption rate p and exhibits a continuous wetting transition at 
a critical growth rate qdp)- The wetting transition is related to the unpinning process of an interface from a substrate 
and may be described by the KPZ equation The model has then been generalized to include a short-range 
interaction between the interface and the substrate . This was done by introducing a modified growth rate go at the 
substrate level. This results in a contact interaction between the interface and the substrate, which is attractive for 
qo < q and repulsive for qg > q. It was found that sufficiently strong attractive interaction modifies the nature of the 
wetting transition, making it first order. In addition, it has been demonstrated that there exists an extended region 
in the phase diagram, where the pinned and the moving phases coexist in the sense that the transition time from the 
pinned to the moving phase grows exponentially with the system size so that the two phases become stable in the 
thermodynamic limit. It should be emphasized that this kind of phase coexistence, which has also been observed in 
the past in other models [ol lid ITU Il2||. can only occur in nonequilibrium systems. 

Some of these results have been confirmed by numerical and mean field studies of KPZ type models [Ulil. In 
particular, by adding a repulsive interaction between the interface and the substrate, a crossover from a continuous 
to a first order wetting transition was found. Non-equilibrium wetting phenomena have also been studied recently in 
models of growing magnetic domains . 

In this paper we present a detailed account of the phase diagram of the SOS model in 1-1-1 dimensions. In order to 
get some indication on the behavior of the model in higher dimensions we introduce a mean field approximation for 
the model, and study the resulting phase diagram and the dynamical behavior of the interface. It is found that the 
mean field approximation reproduces the main qualitative features of the phase diagram obtained in 1-1-1 dimensions. 
In addition, a new feature is found in the pinned phase. By studying the height distribution of the pinned interface, 
two distinct types of behavior are discovered upon varying the dynamical parameters of the model: in one regime the 
distribution decays superexponentially with the height, while in the other the decay is exponential. The two regimes 
are separated by a prewetting transition line. Relevance of this phase diagram to wetting phenomena in dimensions 
higher than 1 -I- 1 is considered. 

The paper is organized as follows. In Sec^Jwe recall the definition of the SOS model and summarize its properties. 
The special case p = 1, where the model is exactly soluble, is discussed in detail in Sec. IIIII The mean field 
approximation, which is the main focus of the present work, is presented and analyzed in Sec. IIVI The paper ends 
with concluding remarks in Sec. 



FIG. 1: Deposition and evaporation processes in the bulk. At the hard-core wall at zero height (not shown here) evaporation is forbidden 
and the deposition rate q is replaced by a modified deposition rate go in order to take the interaction between substrate and surface layer 
into account. 



II. DEFINITION AND PROPERTIES OF THE MODEL 

The SOS model defined in Q is probably the simplest model which follows the spirit of Eq. It is defined on a 
one-dimensional lattice with N sites and periodic boundary conditions, where each site i is associated with an integer 
variable hi describing the local height of the interface. The repulsive part of the potential V{h) is implemented as 
a hard-core wall at zero height, restricting the heights hi to be non-negative. Thus hi — 0,1,2, . Moreover, an 
effective surface tension is introduced by imposing the restricted solid-on-solid (RSOS) condition 

\K-h,±i\<l. (6) 

The model evolves random-sequentially by choosing a random site and carrying out one of the following processes 
(see Fig.!!}: 

- Deposition of an atom on the substrate with rate qq: 

h, = 0^h,^l (7) 

- Deposition of an atom on top of already deposited islands with rate q: 

h, -^h. + l if /ij > 1 (8) 

- Evaporation of an atom at the edge of a terrace with rate r: 

hi -> min(/i,_i, hi,hi+i) (9) 

- Evaporation of an atom in the middle of a plateau with rate p: 

hi ^ hi ~l if hi^i =hi= hi+i > (10) 

If the selected process would violate the restrictions hi > or \hi — hi±i \ < 1, the attempted move is abandoned and 
a new site i is selected. Each attempted update corresponds to a time increment At = 1/N. Since one of the four 
rates can be chosen freely by rescaling time, we set r = 1. Thus the model is controlled by three parameters, namely, 
a growth rate q, a desorption rate p, and a special growth rate qo at the substrate which accounts for an additional 
short-range interaction between the substrate and the wetting layer. 

The model can be easily generalized to higher dimensions. Note that in this case it is possible to introduce different 
evaporation rates for various types of edge sites, e.g. linear edges and corner sites. For simplicity we will assume that 
all rates for evaporation at edges are equal to 1 . In this case the model can simply be generalized to higher dimensions 
by including all nearest neighbors in Eqs. © and ifTUI). 



A. Properties for go = 9 

Let us first consider the case without interactions between substrate and wetting layer, i.e., go — Q- In this case 
the presence of a hard-core wall at zero height leads to a continuous phase transition between a bound and a moving 
phase (see Fig. The phase transition line qdp) is determined by a vanishing propagation velocity of a freely 



4 




FIG. 2: Phase diagram of the wetting model for go = <?■ The second-order wetting transition is represented as a sohd line. The dashed 
line indicates where the coefficient A of the nonlinear term in the KPZ equation effectively vanishes. For p = 1 the dynamical rules obey 
detailed balance (see text). 
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FIG. 3: Example of a closed cycle of deposition and evaporation processes, showing that detailed balance can only be established if p = 1. 



evolving interface. For q < Qc the interface moves downward until it fluctuates close to the wall while for q > qc the 
propagation velocity is positive and the interface detaches from the wall. Obviously this transition takes place even in 
finite systems with a critical threshold depending on N. We note that for p = 0, where evaporation from the middle 
of plateaus is forbidden, the interface velocity in the bulk cannot be negative so that in this case the transition relies 
on a different mechanism 

One can easily verify that the dynamical rules (f7jl- (|10|l do not generally satisfy detailed balance. To this end we 
consider the closed cycle of deposition and evaporation processes, as shown in Fig. 13 Let the statistical weights of 
these configurations in the steady state be Pa, Pb, Pc, and Pq. Obviously detailed balance implies 

2 

q^PA = qPs = Pc^ -Pd - —Pa ■ (11) 

p p 

These equations can only be satisfied if p = 1 (as already anticipated, r is assumed to be equal to 1 without any loss 
of generality). This special case will be discussed in more detail in Sec. IIIII 

The two parameters p and q can be used to control dV/dh and A in the KPZ equation In the case of detailed 
balance, where a bound interface is thermally equilibrated, the coefficient A is expected to vanish. This can be verified 
by comparing the propagation velocities of a horizontal and an artificially tilted interface far away from the wall. The 
line, where both velocities coincide, is shown as a dashed line in Fig. |21 As expected, it crosses the phase transition 
line exactly at the point p = q = 1, where detailed balance is satisfied. 

Moving away from this point, we can therefore study the crossover from equilibrium to non-equilibrium wetting. 
Numerical simulations suggested that the transitions for p < 1 and p > 1 are associated with different sets of critical 
exponents 01 ■ These findings are in accordance with results obtained by Munoz and Hwa j^, who showed that the 
scaling properties of a KPZ interface interacting with a wall depend on the sign of A. 
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FIG. 4: Phase coexistence region in the phase diagram for small values of qo. 



B. Properties for qo < q 



In most experimental applications the wetting layer interacts with the substrate, giving rise to an additional short- 
range force at the bottom layer which may be either attractive or repulsive. In the present model such a force can 
be taken into account by introducing a different growth rate go for deposition at zero height. The influence of this 
parameter was studied in detail in Ref. • Since go does not influence the propagation velocity of the interface far 
away from the wall, the location of the transition line, along which the moving phase is stable, remains unchanged. 
However, if go is smaller than a certain threshold gg, the attractive interaction is so strong as to stabilize the bound 
phase even when the free interface would grow. In other words, for very low values of go, there exists an extended 
region in the phase diagram where the bound and the moving phase coexist in the sense that the transition time from 
the bound to the moving phase grows exponentially with the systems size. The upper boundary of the coexistence 
region depends on go, as shown in Fig.^ 

A thermodynamically stable coexistence of the bound and the moving phase requires a robust mechanism which 
eliminates large protruding islands in the bound phase. In the present model this mechanism works as follows. Once 
an island has been formed by fluctuations, the detached part of the interface quickly grows since q > qc- In the 
phase coexistence region, where the coefficient A in the KPZ equation is negative, the island will grow until the 
slope at the edges exceeds a critical value, where the growth is compensated by the nonlinear term. Afterwards, the 
pyramidial island shrinks linearly with time until it is eliminated. Therefore, phase coexistence can only occur under 
non-equilibrium conditions in those regions of the phase diagram where A is negative. Very recently, the tricritical 
point and the critical behavior at the upper boundary of the phase coexistence region has been investigated in a 
discretized KPZ equation with a potential p^ . 



We first consider the special case p — 1, where detailed balance is satisfied. In this case the stationary proba- 
bility distribution of the bound interface configurations {hi, . . . , h^} is given by the canonical ensemble expression 
corresponding to an energy functional Ti. 



III. EXACTLY SOLUBLE CASE: p = 1 



A. Detailed balance and transfer matrix approach 




(12) 
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FIG. 5: Schematic drawing of V{h) with a potential well at zero height. 



The partition sum 



Zn= exp[-7^(/ii,...,/ijv)] (13) 



runs over all interface configurations obeying the RSOS constraint The energy Ti. is given by 



N 



where V{h) is a potential of the form 

{oo if ft, < 

-Hq/qo) if h^O (15) 

-h ln{q) if ft > 

as sketched in Fig. [3 Therefore, Eq. H12|l may be rewritten as 

P(ftl, . . . , ftAT) = Z^^q^^^^ (q/qo)^^'-- ''-"^ ■ (16) 

Note that in this expression \n{q) plays the role of the chemical potential difference A^/ksT between the wetting 
layer and the gas phase. Obviously the transition takes place at gc = 1- 

To verify the validity of this probability distribution, it is sufficient to demonstrate that the dynamical rules obey 
detailed balance with respect to it. In fact, the deposition process (jH)) reduces the probability P by a factor of q 
while the evaporation processes Q and Hll)|l . which both take place with the same rate, increase P by a factor of 1/g. 
Consequently, the probability currents between pairs of configurations compensate each other so that detailed balance 
is satisfied, proving the validity of the equilibrium ensemble (I16II . Similarly one can show that detailed balance is also 
satisfied at the bottom layer. We note that these considerations are only valid in the bound phase q < 1, where the 
probability distribution is stationary. Once the system enters the moving phase, the process is out of equilibrium. 

The canonical ensemble can be used to compute the density profile of a bound interface in the case of detailed 
balance. To this end we use a transfer matrix formalism introduced in [I7lll8j |. writing the Boltzmann factor exp(— 
in Eqs. ifT^ and lfTB|l as a product 



Pihu...,hN) = ^l\n^M+i, (17) 

where Zn = Tr(T^) and 



l-h,h' 



(18) 

otherwise. 
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The transfer matrix T is infinite-dimensional, acts in spatial direction, and yields the contribution to the Boltzmann 
factor between adjacent sites with the heights h and h' . Because of the RSOS condition ||HJl, it has a tridiagonal 
structure and reads 



T - 



\ 



^3/2 ^2 ^5/2 



^5/2 ^3 ^7/2 



V 



(19) 
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Using the transfer matrix formalism, the stationary density p{h) of sites at height h can be expressed as 

p{h) = Z-^\h\T''\h), (20) 

where {{hW denote canonical basis vectors in height space. For iV — > cxd this expression is governed by the 

largest eigenvalue A of the transfer matrix 



T|0)-A|0), 



(21) 



where \(j)) denotes the corresponding eigenvector. Thus, in an infinite system, the stationary densities may also be 
written as 



p{h) = 



(22) 



Note that the numerator in this expression is quadratic in |0), just as in a quantum- mechanical problem. 



B. The case qo — q 



In order to understand how the interface detaches from the wall, it is useful to study the scaling behavior of 
the density of sites at the bottom layer close to the transition point. Let us first consider the case go — q, where 
interactions between substrate and bottom layer are absent. For h > the eigenvalue problem reads 

g"-'/'0/.-i + + 'z"+'/'0/.+i = , (23) 

where A is the largest eigenvalue of T and (j)h are the components of the corresponding eigenvector |0) representing 
the stationary state. Close to criticality we can carry out the continuum limit (j>h <t>{h)^ replacing discrete heights h 
by real- valued heights h. In this limit the above eigenvalue problem turns into a differential equation which, to leading 
order in e = 1 — g > 0, is given by 

(^ + (3-A)-3e/i)(/.(M = 0. (24) 
This differential equation is solved by an Airy function 

(3e)^/3 )■ (25) 

Since (t>(h) has to vanish for h <0, we obtain the eigenvalue A = 3 so that ^(/i) — Ai([3e]^/'^/i). Therefore, the heights, 
in particular the average height and the interface width, scale as 

{h)^wr^e-^^\ (26) 

Next, we determine the density of exposed sites at the substrate p(0). In the continuum limit of Eq. ()22|l (010) is 
proportional to 0'(O), hence 

p(O) = ■l[0'(O)]^ (27) 
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FIG. 6: 1+1-dimensional interface with island sizes ^i,^2i^3- 



where Af — dh (p{h) is a normalization factor. Since 0'(O) - and A/" - £-1/3 one obtains a linear scaling law 

p{0)^e. (28) 

Thus the density of exposed sites at the bottom layer scales linearly with the distance from criticality, proving that 
the wetting transition at p = qc = 1 is continuous. 

Imposing fixed boundary conditions /ii = /iat = 0, it is also possible to study finite-size scaling at the critical point 
Qo — 1 = ^- At criticality the transfer matrix has a simple structure and can be thought of as generating a simple 
random walk near a wall so that the mean height and the bottom layer density scale as 

{h) - , p(0) - (29) 

The transfer matrix does not provide any information regarding dynamical properties. Numerical simulations (whose 
details are not shown here) suggest that an initially flat interface at zero height roughens with time (for t <C N'^) as 

{h) - , p(0) - (3Q^ 

Assuming standard power law scaling, we can combine Eqs. if^ . (j^ . and (|Sn|l in the scaling forms 

{h{e,N,t))=N'/'f(t/N^,eN'/') 

(31) 

where / and g are scaling functions with an appropriate asymptotic behavior. As expected these scaling forms are 
consistent with the critical exponents z = 2, a = 1/2 of the Edwards- Wilkinson universality class 

Another interesting aspect is the stationary distribution P{£) of island sizes £ in the bound phase (see Fig.|Bl). In the 
case of detailed balance this distribution can be computed exactly. To this end we introduce the projection operator 
P — 1 — |0)(0| which projects onto states with nonzero height. Moreover, let Q = FTP be a transfer matrix describing 
an interface that does not touch the bottom layer. Obviously the distribution P{£), which may be interpreted as a 
first-return probability of the interface to the bottom layer, is given by 

_ (o|rg^'^r|o) ^ (o|rg^-^r|o) 
^ ' " (OIT^IO) - A^ ' ^'^ ' 

Note that the leftmost column and the topmost row of the restricted transfer matrix Q are zero, while all other matrix 
elements are the same as in Eq. ifT^ . Because of this simple structure, the spectrum of Q is just the spectrum of T 
multiplied by q combined with a zero mode so that the largest eigenvalue of Q, which dominates the matrix product 
in Eq. (|32|l , is qA. In the limit £ — > oo we therefore obtain an exponential distribution of the form 



Therefore, the average island size scales as 



P{£)-q'. (33) 



(34) 

mq e 



and diverges at the transition. Since p(0) — l/£, this result is in agreement with Eq. (|28|l . 

In order to determine the stationary interface profile p{h) in the limit h oo, let us go back to Eq. (|23|l . By 
assuming that the second and the third term on the left hand side can be neglected, one finds that 

0, ~ i g^-i/^ , (35) 
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FIG. 7: Numerically determined stationary interface profile of the full model in 1+1 (lower data points) and 2+1 dimensions (upper data 
points) for different values of p slightly below the critical threshold. The bold lines indicate the slopes 2 and 3, respectively. 



with the corresponding asymptotic solution 

Therefore, in the case of detailed balance, the profile of a bound interface decays as a Gaussian for large values of h. 
Even for p 7^ 1, where detailed balance is violated, numerical simulations (see Fig.[7)l suggest that in 1+1 dimensions 
the entire bound phase is characterized by Gaussian density profiles. 



C. First-order phase transition for small values of go 



Let us now consider the influence of an attractive force between substrate and wetting layer by taking go < 1- 
Obviously, the interface is bound to the wall for q < 1, while for q > 1 it will tunnel through the potential barrier. 
This means that the critical point q^ — X remains unchanged. However, if go is decreased below a certain threshold gg, 
the attraction is sufficiently strong such that the transition becomes first order. 

To demonstrate the crossover to a discontinuous transition in the case of detailed balance, we look for a localized 
pinned interface solution at the transition point g = 1. Assuming an exponential interface profile 

(i>h for /i > 1 (37) 

with some z < 1, the eigenvalue problem then reduces to three independent equations 

-1 1/2 

% 00 + 90 ^ = 

q^'/^ct>o + z + z' = Az (38) 

z"^ + 1 + z = A 

which have the (unnormalized non-negative) solution 



, 1/2 + - 3ql 1 z + 1 



2(1 -go) 2' go 

Consequently, the stationary density of exposed sites at the bottom layer is 



(n\ - - go _ 1 + go - 6gg + ^1 + 2go - 3g§ 

Eh=o K 90 + 2 + 4go - 6g^ 



while the densities at higher levels are given by 



p{h) P ^ . (41) 
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FIG. 8: Phase diagram for p = 1 in the (g, qo)-pla.ne. 



It turns out that the bottom layer density p(0) is positive for qq < 2/3 and vanishes at qo — 2/3. For qq > 2/3, 
however, one has z > 1, so that the exponential ansatz (j)fi = z'^ is no longer physically meaningful. Hence, in the case 
of detailed balance, the transition becomes first order at the tricritical point 



P = gc = 1 , 



go- 3 



(42) 



For go < §j the potential well is deep enough to bind the critical interface to the wall, leading to an exponentially 
decaying interface profile. For go > |, such a localized solution does not exist and the transition becomes continuous. 



D. Scaling properties near the tricritical point 



The phase diagram for p = 1 is shown in Fig. |H1 Using the transfer matrix approach, we can show that the critical 
behavior along the second order transition line is always the same as the one discussed in Sec. IIII Bl However, in the 
vicinity of the tricritical point the scaling properties are different. Moreover, they depend on the direction from which 
the tricritical point is approached. 

Let us first consider the case g = 1, approaching the tricritical point horizontally from the left along the first-order 
phase transition line. Using the expressions (|40() and H41I) one can compute the interface height 



h = ^MM = - 

/i=0 ^ 

and the squared interface width 

oo 



2go 



go - 6g^ + + 2go - 3g^ 
2(1 + 4go - 3g3 + (g^ - 3go - 1) ^1 + 2go - 3g, 



(43) 



h=0 



l + 3go-3Vl + 2go-3g^ 



(44) 



Approaching the tricritical point from the left by increasing go, these quantities scale to lowest order in (5 = gp — go 

as 



p(0) = 4(5 



(45) 



On the other hand, if the tricritical point is approached vertically keeping go = go = 2/3 fixed, a numerical diago- 
nalization of the transfer matrix suggests that the asymptotic interface profile crosses over from an exponential to a 
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FIG. 9: Decay of the density of sites at zero height at the tricritical point q = p = 1, qo = 2/3. The slope of the curve tends to —0.24(1), 
leading to the conjecture that p(0) ~ t"^/'*. 



Gaussian decay. In this case the height, the width, and the bottom layer density are found to scale as 

where e = I — q. Moreover, by keeping the boundary sites fixed at zero height, we can study finite-size scaling at the 
tricritical point. Evaluating products of the transfer matrix numerically we find that 

h ^ 

(47) 

Finally, numerical Monte Carlo simulations at the tricritical point (see Fig. O suggest that an initially flat interface 
roughens in such a way that 

h ^ 

Assuming standard power law scaling all these results can be combined in the scaling forms 

h{e, S, N, t) = N^'^F{t/N^,eN^'^, 5N^'^^ 
po{e, S, N, t) = N-^'^G{t/N^,eN^'^, 5N^I^ 



(49) 



where F and G are scaling functions with an appropriate asymptotic behavior. Again these scaling forms are consistent 
with the critical exponents z = 2, a = 1/2 of the Edwards- Wilkinson universality class. 



IV. MEAN FIELD THEORY FOR NON-EQUILIBRIUM WETTING 



Mean field theories describe a system in an approximate way by ignoring spatial correlations. For a given model 
there are many possible types of mean field theories, depending on the microscopic level one is trying to describe. 
Previous mean field approaches to non-equilibrium wetting considered the average height as the dynamical variable 
and studied the mean field approximation of its dynamics. For example, in a study by Giada and Marsili the KPZ 
equation Q was mapped by a Hopf-Cole transformation to a Langevin equation with multiplicative noise, discretizing 
space and replacing nearest-neighbor interactions by global couplings. Using a Morse potential with a potential well 
at zero height they were able to reproduce second- and first order transitions as well as phase coexistence. More 
recently, Santos et al. 0| extended these studies, suggesting that the mean field phase diagram does not change if 
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fluctuations are taken into account. Moreover, they identified a narrow domain close to the borderline between phase 
coexistence region and wet phase, where the system exhibits spatio-temporal intermittency. 

In this Section we construct mean field equations describing the temporal evolution of the height distribution ■0„ 
for height h = n. This approach is expected to yield more detailed information on the structure of the interface than 
previous mean field theories. 



A. Mean field equations 



In order to construct the mean field equations for the height distribution, let us first consider the case of a freely 
evolving interface without a wall. The rate equations describing the temporal change of V'n consist of several terms 
corresponding to different processes. Let us, for example, consider the probability that an attempted update leads 
to the desorption of an atom from the interior of a plateau at level n. This probability is proportional to ipn, which 
is the probability to find a randomly selected site at height n. Moreover, it depends on the heights of the nearest 
neighbors, which are restricted to take the values {n — 1, n, n + 1}. For simplicity we assume that each site has two 
nearest neighbor sites. Clearly one can generalize it to the case where each site has 2{d— 1) nearest neighbors. In this 
case one may have several types of edge of corner sites, requiring a larger number of growth rate parameters. To avoid 
this complication, we restrict ourselves to the case of two nearest neighbors. Thus, neglecting spatial correlations, the 
probability to find the two nearest neighbor sites at the heights Z,m £ {n — 1, n, n + 1} is assumed to be proportional 
to ipi'4'm divided by (V'n-i + V'n + V'n+i)^- For example, for evaporation from a plateau we have n = I ^ m and thus 
the contribution of the process H10() to the dynamical equation of ipn is 

(50) 



{lpn-l + V'n + Ipn+l)"^ ' 

Similarly the loss of probability due to evaporation at the edges © is given by 
while the depositon process (O leads to a loss term of the form 



(51) 



^ (V'n-l +V'n+l)^ 

In addition, there are corresponding gain contributions at the neighboring levels such that the total probability is 
conserved. Collecting all terms, the mean field equations can be written as 

— TT — ^ J (53) 

at 

where 

{tpn + i^n+l + ^71+2^ {tpn-l + tpn + '4'n+iy 

The hard-core wall at the bottom layer can be taken into account by formally setting ip-i = A-i = and replacing 
q with qq in the expression for Aq, which then can be written as 



B. Mean field phase diagram 



The main result of the calculations, which will be presented in detail below, is the mean-field phase diagram shown 
in Fief. 1101 for the case go = As in the 1 + 1 dimensional model, there is a continuous phase transition from a bound 
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FIG. 10: Mean field phase diagram for go = 9- The second-order phase transition line is represented as a solid line. The bound phase 
consists of two parts, where the interface profile ipn either decays expontially (E) or superexponentially (SE) for large n. 

to a moving phase. In calculating the mean field height distribution in the moving phase it is found that the interface 
is localized around its average height at any given time, representing a smooth growing interface. This is expected in 
high dimensions, for which mean field represents a reasonable approximation. Clearly, in d = 1 + 1 dimensions this is 
not the case as discussed before. 

In the bound region, two types of phases were found. In the larger part of the p, g-plane (denoted as SE in Fig. llU|) 
the height profile decays superexponentially at large height. In particular it is found that for large n the profile takes 
the form 

i>n ~ g"" exp(-2"-") , (56) 

where a is a constant. In addition, another region is found (denoted as E in Fig. I10|l . where the height distribution 
decays exponentially 

V'n ~ e-"" , (57) 

where a > is a constant. 

The superexponential behavior may be understood by considering the equilibrium case (i.e. p ~ 1) in which a 
d — 1-dimensional manifold is attracted by a gravitational force to a hard wall. The energy of a fluctuation reaching 
a height n scales as n'^, leading to the height distribution 

^ e-'^"" (58) 

where /? > is a constant. This behavior is seems to be valid in an extended non-equilibrium region of the phase 
diagram, as shown in Fig. \7\ where numerical simulations in d = 2 + 1 dimensions are presented. 

Numerical attempts to find signatures of a possible exponential phase in 2 and 3 dimensions failed since it is very 
difficult to obtain a reliable statistics in the tail of the height distribution, especially in higher dimensions, where 
finite size effects become increasingly relevant. The conjecture, that the exponential phase might be related to the 
roughening transition of KPZ interfaces in d > 2 could not be substantiated by numerical simulations. 



C. Mean field equations 



In the stationary state one has Aq = and An — A„_i = so that An vanishes for all n > 0. Therefore, the 
stationary mean field equations read 

= p^l + 2^-0^? + ^0^1 - qo Mi^Q + V^i + ^2)^ , (59) 

= p-0n-H + 2-0nV'n-H + V^nV^"+l - g-0n(V^« + V'n-FO^f ^" ^"t^ ^ ^"^^ ) , (60) 
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FIG. 11: Classification of fixed points depending on p and q. In region I (II) there exists only one real fixed point x* < 1 {x* > 1), 
while in region III there are three real fixed points. Regions I and II are separated by the straight dotted line q = {p + 3)/4. The phase 
transition line (solid line) crosses from region I to region II and back into region I, as illustrated in the upper inset. The lower inset zooms 
the area where the both lines enter wedge-shaped region III. 



where n=l,2,...,oo. In order to solve this equation by iteration, it is convenient to consider quotients of successive 
densities 



1pn-l 

In terms of these variables, the stationary mean field equations take the form 
= pxl + 2x1 + Xi - qo{l + a;i + xiX2f , 

n 3 , o 2 , n , ^2 2('^+Xn+l+Xn+lXn+2\'^ 

= px„+i + 2x„+i + Xn+1 - q{l + Xn+l) xA—— ■ . 

\ l+Xn+ XnXn+1 ' 

The bulk equation (|63|) can be interpreted as an iterative map (a;„, a;„-|-i) (a^n+i, 2^n+2), where 



Xn^2 — f{XmXn+l) 



-1 



1 



1+Xn+ XnXn+l / 1 + 2x„+i + px^^^ 



(l + a;„+i) y qxn+i 
for n = l,2,...oo. In addition the initial condition for xi, X2 is given by the bottom layer equation 

1 ll + xi{2+pxi) 



X2 



-1 



Xl 



qoxi 



(61) 

(62) 
(63) 

(64) 

(65) 



D. Stationary solutions 



In order to evaluate the stationary height distribution for given p and g, one has to look for a point on the line 
of possible initial conditions (|65|l and iterate the map H64|l such that it reaches a real fixed point with x* < \. This 
trajectory then corresponds to a physical height distribution. 

To proceed, we first analyze the fixed points of the map. The fixed point equation 

X {qx^ + {2q - p)x^ + {q - 2)x - =0 (66) 
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has four solutions. Two of the solutions, a;g = and xl are real in the entire p, q-plane, while other two solutions 
X2,x^ may either be both real or complex conjugate to each other. We denote the region in the p, g-plane, where 
X2 J are real, by region III (see Fig. 111(1 . We further divide the complementary region to III into two regions, I 
where a;* < 1, and II, where a;* > 1. Since the bracket in Eq. IjtiGI) is equal to —1 at x = and positive for g > in 
the limit of large x the fixed point xl has to be positive. 

Analyzing the map ((64|l with the initial condition H65|l we find that the only fixed points, which correspond to 
physical height distributions, are either Xq or when it is real. Here is the solution of Eq. (|65|l which satisfies 
< 1 < 0:3 in the region III. In particular we find that below the transition line, the height distribution is controlled 
by the fixed point = in region I and by the fixed point X2 in region III. Details of the analysis, which led to this 
result, are given in the Appendix. 

We now study the stationary height profile in the two regions. In region I the map flows to the hyperbolic flxed 
point a;^ = along its stable trajectory. Expanding the map for small values of x to lowest order, this stable manifold 
is described by the nonlinear relation 

Xn+i^qxl, (67) 
in the limit a; — > 0. Along this manifold, the map approaches the fixed point a;Q = super-exponentially as 

x„ ~ iexp(-2"-") (68) 

yielding the height profile in Eq. (|56|l . 

In region III stationary solutions are controlled by the fixed point X2 > (see Appendix). Linearizing the map 
around this fixed point one obtains an exponential behavior 

Xn e-™ , (69) 

where a > is a constant, leading to the exponential height profile H57|l . 

To complete the analysis of the phase diagram, one has to locate the wetting transition line. This is done by 
simulating the dynamical mean field equations H53|l and determining the point from where on the interface detaches. 
This analysis leads to the transition line shown in Fig. 1101 

So far all mean field results were obtained for Qq = q, i.e., without an attractive force between the substrate and 
the wetting layer. Lowering go changes the bottom layer equation (|65l) and thereby the possible starting points of the 
iteration. As shown in the Appendix, the mean field approximation reproduces the phenomenon of phase coexistence. 
However, unlike in the model in 1 + 1 dimensions, it emerges everywhere along the phase transition line as soon as 
go < 9; hence, there is no tricritical point. Moreover, in the limit qo the threshold for the growth rate q, where 
the interface detaches, tends to infinity. Therefore, the region of phase coexistence is not bounded from above as in 
the 1 + 1 dimensional model. 



V. CONCLUSIONS 



In this paper we have given a detailed account of a recently introduced solid-on-solid model for non-equilbrium 
wetting, which is defined in the spirit of a KPZ equation in a potential. Introducing a hard-core wall at zero height, 
the model exhibits a contininuous wetting transition from a bound to a moving phase. The model is controlled by 
two paramters p and q, which effectively determine the asympotic slope of the potential and the coefficient of the 
nonlinear term in the KPZ equation. 

For p = 1 the dynamical rules of the model obey detailed balance. In this case the stationary distribution of a 
bound interface is given by a Boltzmann ensemble, which allows one to derive various quantities exactly. Moving 
away from this line, the model crosses over to a non-equilibrium behavior, that is characterized by different critical 
properties. The model can be generalized further by including an attractive interaction between the substrate and the 
wetting layer. If this force is strong enough it may turn the continuous transition into a discontinuous one. Moreover, 
the bound and the moving phase may coexist in regions where the coeffient of the nonlinear term in the KPZ equation 
is negative. 

In order to assess the behavior of the model in higher dimensions, we have proposed a mean field approximation, 
which is based on rate equations for the densities at different heights. The phase diagram turns out to be surprisingly 
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rich. It turns out that the mean field approximation reproduces the properties of the original model. However, in the 
moving phase the interface remains smooth as it is expected for KPZ-type growth in higher dimensions. 

As a new feature, the mean field approximation predicts the existence of two different regions in the bound phase. In 
one of these regions the interace profile decays superexponentially with increasing height while in the other region an 
exponential decay is observed. This can be explained by classifying the fixed points of an iterative map for quotients 
of the densities. The question to what extent this crossover from superexponential to exponential profiles can be 
observed in the full model is still open. 
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FIG. 12: Stationary solutions in region I: Superstable manifold of the fixed point Xq = for p = 2 and different values of q in regions 
I and II. For q < qc the manifold (bold line) interesects with the with the dashed line of possible initial conditions I65i . representing a 
stationary solution in the bound phase, where the bottom layer density ipo is positive. Approaching qc this intersection point moves to 
infinity and tpg tends to zero. For q > qc the superstable manifold originates in the other real fixed point so that no stationary solution 
exists. Similar graphs are also obtained for other values of p in regions I and II. 



APPENDIX A: STATIONARY SOLUTIONS OF THE MEAN-FIELD EQUATIONS 



1. Fixed points 



Besides Xq = 0, Eq. has three other fixed points which are the roots of the polynomial equation 

qx^ + {2q-p)x^ + {q-2)x-l = 0. (Al) 

As shown in Fig. 1111 the p, g-plane can be divided into three different regions. In regions I and II one fixed point is 
real and two of them are complex conjugate, while in region III all fixed points are real. The boundary of region III 
is characterized by the existence of a two- fold degenerate fixed point. Comparing Eq. (jAlll with a polynomial of the 
form (x — a)'^{x — b) ~ 0, the boundary of region III can be given in a parameter representation by 

q = J , P = 2 — , b = —— , (A2) 

where < a < 1. The boundary has the form of a wedge, shown as a dashed line in Fig. The lower branch for 
< a < 1/2 and the upper branch for 1/2 < a < 1 terminate in the triple point (p, q) — (8, 28), where Eq. IjAljl has 
a three- fold degenerate fixed point x* = 1/2. 

Physically meaningful stationary solutions of the iterative map must start from a point on the line given by Eq. H62|l 
and have to fiow towards a real fixed point with < x* < 1. For this reason we divided the complement of region III 
into two parts, namely, region I with x| > 1 and region II with x^^ < 1. Both regions are separated by a straight line 
q= {p + 3)/4, where xl = 1. 



2. Stationary solutions in regions I and II 



We first show that in regions I and II a physically meaningful stationary solution always flows to the fixed point 
Xq = 0. To this end we show that the other real fixed point x* is either larger than 1 or unstable. 

As shown in the upper inset of Fig. 1111 the phase transition line starts at p = in region I, crosses into region 
II at the point q = p = 1 and then crosses back into region I at the point {p,q) ~ (5.380,2.095). Obviously the 
fixed point x\ can only be physically meaningful between these two crossing points, where x| < 1. However, in this 
interval x^ turns out to be unstable. To demonstrate this point we consider the eigenvalues of the Jacobian of the 
map Xn+2 = fixn,x„+i) in Eq. 



Ai,2 - I 



df 

dXn+ 



± 



- df 

dXn 



+ 1 



OXn 



(A3) 




FIG. 13: Fixed points (bullets) and possible stationary solutions (bold lines) for p = 60 inside for (A) q = 12.0, (B) q = 13.5, and (C) 
q = qc ~ 15.644. The figure is explained in the text. 



Using Eq. (|A1|I the partial derivatives can be expressed as 

df 1 



df 



qx 



2qx*^ + 2{q - l)x*'^ + {3q - 2)x* - 2 



(A4) 



Along the dotted line in Fig. where x* = 

Ai,: 



2qx*^{x* + 1) 

1, the two eigenvalues 
1 



1- -(3 ±79^^249) 



(A5) 



are complex conjugate in the interval between the two crossing points 1 < g < 5.380. As can be verified numerically, 
they are also complex conjugate in a neighborhood of this line, which includes the phase transition line. Since the 
determinant of the Jacobian 



A1A2 - — o 



(A6) 



is larger than 1 for x* < 1, the fixed point x* is found to be unstable. Thus we can conclude that physically meaningful 
stationary solutions in regions I and II are always controlled by the fixed point Xq = 0. 

Since for a;* ^ the two eigenvalues tend to Ai = and A2 — —00, the fixed point Xq = is nonlinear and 
hyperbolic. It has superstable manifold which in the limit of x ^ is given by a;„+i = qx^. This manifold is shown 
in Fig. 1121 for different values of q below, at, and above the critical point. As it can be seen, a stationary solution 
exists if the manifold intersects the line of possible initial conditions H65|). 



3. Stationary solutions in region III 



In the wedge-shaped region III, there are three real fixed points x^ < X2 < x"^. The first one is smaller than 1 and 
unstable, while the second one is smaller than 1 and hyperbolic. The properties of Xg depend on p and q. Below the 
dotted line in Fig. ^2 it is larger than 1 and stable, while it is smaller than 1 and unstable above. 

Fig. El illustrates typical situations at p = 60 for various values of q, crossing the region III from bottom to top. 
Below the wedge in region I (panel A), the stationary solution can be constructed in the same way as before. Entering 
region III from below (panel B), two new real fixed points x| < X2 < 1 emerge, which are fully unstable and hyperbolic, 
respectively. The superstable manifold of the fixed point Xq = now originates in while it is the linearly stable 
manifold of the hyperbolic fixed point Xj which intersects the dashed line of possible initial conditions (|65|l . As before, 
the intersection point moves continuously to infinity as q approaches the critical point (panel C). Thus the unbinding 
transition manifests itself in the same way as in the regions I and II, the only difference being that the physically 
relevant stable manifold is now controlled by the hyperbolic fixed point X2 > instead of Xq = 0. The linear stability 
of X2 is responsible for the purely exponential profile observed in region E. 
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FIG. 14: Phase coexistence in the mean field approximation for p=0.8 (see text). 



4. Phase coexistence in regions I and II 



Lowering qq changes the hne of possible initial conditions l|t)5|l . while the stable manifold of the fixed point Xq = 
remains the same. The typical situation for p = 0.8 is shown in Fig. 1141 The left panel shows the stable manifold (solid 
line) and the bottom layer equation (dashed line) without attractive force at the critical point go = 9 = 9c — 0.943. 
Both curves approach each other smoothly and intersect at infinity so that the transition is continuous. The right 
panel shows the same situation in the presence of an attractive force for qq = 0.5. Accordingly, the effect of lowering 
qa is to move the dashed line of possible intial conditions upward such that it intersects the critical stable manifold 
at a certain finite point. This means that the bottom layer density ipo is finite at the critical point, making the 
transition first order. Moreover, a stationary solution still exists even if q is increased, proving the possibility of phase 
coexistence in the mean field equations. Increasing q beyond a certain threshold, a stationary solution no longer 
exists. This defines the upper boundary of the phase coexistence region in the mean field phase diagram. 
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